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Light scalar hadrons can be understood as dynamically generated resonances. These arise as 
‘companion poles’ in the propagators of quark-antiquark seed states when accounting for meson- 
loop contributions to the self-energies of the latter. Along this line, we extend previous calculations 
of Tornqvist and Roos and of Boglione and Pennington, where the resonance do (980) appears as 
companion pole in the propagator of do (1450) which is predominantly a quark-antiquark state. We 
also construct an effective Lagrangian where do (1450) couples to pseudoscalar mesons with both 
non-derivative and derivative interactions. Computing the one-loop self-energy, we demonstrate 
that the propagator has two poles: a companion pole corresponding to do (980) and a pole of the 
seed state do (1450). The positions of these poles are in quantitative agreement with experimental 
data. 
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I. INTRODUCTION 

The majority of mesons can be understood as being 
predominantly quark-antiquark states [T|. Yet, various 
unconventional mesonic states such as glueballs, hybrids, 
and four-quark states are expected [2J. One particular 
type of four-quark meson is that of a ‘dynamically gen¬ 
erated’ state Although there are recent works on 

the compositeness of (dynamically generated) resonances 
- see e.g. Refs. mama, where the former also provides 
a method to quantify the weight of the two-body com¬ 
ponents of a resonance - a generally accepted definition 
of dynamical generation does not exist inj- An interest¬ 
ing version of this idea was put forward in Refs. huhh. 
Consider, for instance, a single seed state, e.g. a quark- 
antiquark meson with certain quantum numbers. This 
state interacts with other mesons, giving rise to loop con¬ 
tributions in the corresponding self-energy. These shift 
the pole of the seed state and, moreover, a(t least one) 
new pole may appear. The latter, also denoted as a com¬ 
panion pole, corresponds to a dynamically generated res¬ 
onance. As a consequence, two resonances have emerged 
from a single seed state. In this work, we aim to review 
some previous works and deepen the understanding of 
the resonance a o (980) as a dynamically generated state. 

There is a growing consensus that the scalar reso¬ 
nances /o(1370), /o(1500), i\o(1430), and ao(1450) are 
predominantly quark-antiquark states; see for example 
Refs. [U 0Q31 ITS1 - [7T| [recent studies [721H3] agree in in¬ 
terpreting the resonance /o(1710) as predominantly glu- 
onic]. Then, the light scalar states /o(500), /o(980), 
A'g(800), and ao(980) are (most likely) predominantly 
four quark-states [see e.g. Refs. [SHU [II na MM and 
refs, therein]. As we shall show in the following for the 
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heavy scalar-isovector seed state a o (1450), the coupling 
of this state to nr], I\K , and nrf dynamically generates 
the light ao(980) as a particular type of four-quark me¬ 
son. 


Tornqvist and Roos El (in the following denoted as 
TR) and later Boglione and Pennington |T5J jT3] (denoted 
as BP) studied the mechanism of dynamical genera¬ 
tion through meson-loop contributions to the self-energy. 
Here, we extend their studies (Sec. 0 and compare nu¬ 
merical results for the poles of the propagator to the lat¬ 
est experimental data pQ. It turns out that (depending 
on the assignment of the poles to physical resonances) 
the widths of both the seed state ao(1450) and the dy¬ 
namically generated state ao(980) are by a factor of 2 
larger than the experimental values. Moreover, the mass 
of a o (1450) is too large (by 100 MeV in TR and by 400 
MeV in BP). It thus seems that, while qualitatively fea¬ 
sible, the dynamical generation of resonances as compan¬ 
ion poles in the propagator does not yield results that are 
in quantitative agreement with experimental data. 


In this work, we show that this is actually not true and 
that the mechanism of dynamical generation produces re¬ 
sults which are in quantitative agreement with the data. 
To this end, we introduce a Lagrangian inspired by the re¬ 
cently developed extended Linear Sigma Model (eLSM) 
[20l 21 . Here, the mesons interact via derivative and 
non-derivative couplings (Sec. Ill). In our case, the La¬ 
grangian contains a single scalar-isovector ao seed state 
which corresponds to the resonance ao(1450). A careful 
analysis of the pole structure of the corresponding prop¬ 
agator shows that it is indeed possible to obtain a narrow 
resonance with mass around 1 GeV, the pole coordinates 
of which fit quite well with those of the physical ao(980) 
resonance, and simultaneously obtain a pole for the seed 
state in agreement with that for the ao(1450) [T]. Fi¬ 
nally, we also investigate the use of Feynman rules in the 
context of quantum field theories with derivative interac¬ 
tions, and demonstrate that for a particular form of the 
Lagrangian there may be a discrepancy between ordinary 
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Feynman rules and dispersion relations (see App. @- 

Our units are H = c = 1 . The metric tensor is r]^ v = 
diag(+, 

II. DYNAMICAL GENERATION 
A. Approach of TR and BP 

Following earlier work jlOj . Tornqvist et al. studied the 
scalar sector in a unitarized quark model by including 
meson-loop contributions m- They showed that meson- 
loop effects may serve to explain the existence of the light 
scalar states. 

The following two points are relevant in the mechanism 
of dynamical generation, irrespective of the quantum 
numbers of the hadronic resonance considered: (i) The 
propagator of a quark-antiquark seed state gets dressed 
by meson-loop contribution to the self-energy. These con¬ 
tributions shift the mass of the state and change the form 
of its spectral function. When increasing the coupling, 
the corresponding pole moves away from the real axis 
and follows a certain trajectory in the complex plane. 
The mass and the width of the resonance are determined 
by the position of the complex pole of the dressed prop¬ 
agator on the appropriate Riemann sheet - a procedure 
first proposed by Peierls a long time ago [30]. (ii) If the 
interaction exceeds a critical value, a(t least one) com¬ 
panion pole can appear in the complex plane. If this pole 
is sufficiently close to the real axis, it can manifest itself 
in the spectral function as an additional resonance with 
the same quantum numbers as the seed state |2I ED EH- 
Since the coupling of scalars to pseudoscalars is large, the 
scalar sector is particularly affected by such distortions 
of the spectral function. 

We now recapitulate the seminal works TR m and BP 
m, where the latter uses the same model as the former 
but with a slightly different set of parameters. The main 
goal is the determination of the inverse propagator of a 
resonance after applying a Dyson resummation of loop 
contributions to the self-energy: 

A -1 (s) = s — ml — n(s) , (1) 

where s is the first Mandelstam variable, mo is the bare 
mass of the seed state, and II(s) = JV IR (s) is the self¬ 
energy [32!. Here, the sum runs over the loops emerging 
from the coupling of the resonance to various mesons. 
The imaginary part of n^s) corresponds to the partial 
decay width of the resonance into mesons in channel i. 
The real part of n(s) on the real axis is related to the 
imaginary part by the dispersion relation 

Ren(s) = -Ids' ~ ImII( / S, ' > . (2) 

7T J S — S' 

TR and BP now assume a simple model for the imaginary 
part of n, (s), see Refs. ESOE3] for details: 

IrnffR) = — sa,i)F^(s)Q(s — s t h,i) ■ (3) 

Vs 


In the scalar-isovector sector the Adler zeros SA,i are set 
to zero for simplicity mm- The form factor is chosen 
to be a simple exponential, 

Fi(s) = exp[— kl{s)/(2kl)\ , (4) 

where is a cutoff parameter, and ki(s) is the absolute 
value of the three-momentum of the decay particles in 
the rest frame of the resonance, 

h (s) = 2^/! \J s2 + ( TO ?i - m i 2 ) 2 ~ 2 Ofi + «4)s ■ (5) 

Here, mu, are the masses of the decay particles, i.e., 
in our case the pseudoscalar mesons [34] . The function 
Fi(s) guarantees that the imaginary part of n(s) vanishes 
sufficiently fast for s — > oo (the inverse cutoff fco corre¬ 
sponds to the non-vanishing size of a typical hadron). 
The step function in Eq. (3j) ensures that the decay chan¬ 
nel i contributes only when the squared energy of the 
resonance exceeds the threshold value Sth,i • Finally, the 
coupling constants gi are related by 5'17(3)-fiavor sym¬ 
metry. 

Note that one may also define the so-called Breit- 
Wigner mass of a resonance as the real-valued root of 
the real part of the inverse propagator, ReA _1 (s) = 0. 
These roots can be found by identifying the intersections 
of the so-called ‘running mass’ 

m 2 (s) = TOq + Ren(s) (6) 

with the straight line /(s) = s, where s is purely real. 
This definition of the mass of the resonance is also used 
in TR and BP. However, the Breit-Wigner mass does not 
necessarily correspond to a pole in the complex energy 
plane or to a peak in the spectral function. 

For the scalar-isovector sector, the main results of TR 
and BP can be summarized as follows: 

1. TR found a pole on the second Riemann sheet 
with coordinates [3S] m po i e = 1.084 GeV and 
Fpole = 0.270 GeV, which is a companion pole cor¬ 
responding to the resonance ao(980). A reanalysis 
[the second paper quoted in Ref. [IT] ] where the 
complex plane was investigated more carefully re¬ 
vealed another pole with TO po i e = 1.566 GeV and 
Fpoie = 0.578 GeV on the third sheet. This pole is 
indeed the original seed state and describes the res¬ 
onance ao(1450). It was suggested that, although 
the numerical agreement was not yet satisfactory, 
an improved model could in principle be capable 
of describing the whole scalar-isovector sector up 
to 1.6 GeV. TR also reports one (but not more) 
intersection point (s) of the running mass from Eq. 

<©• 

2. BP used the same approach, but did not look for 
poles of the propagator. Instead, they considered 
the Breit-Wigner mass. Compared to TR, also the 
values of the bare mass parameter toq as well as 
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the overall strength of the couplings gi in Eq. ^ 
were changed. BP found two intersection points 
for the running mass from Eq. ([6]), one in the re¬ 
gion around 1 GeV corresponding to ao(980) (like 
TR) and another one at about 1.4 GeV (absent in 
TR). This latter intersection was interpreted as the 
state Oo(1450). Note that, although BP did not in¬ 
vestigate the poles of the propagator, a pole and an 
intersection were reported in an earlier work [12] . 

Apparently, the situation is somewhat inconclusive re¬ 
garding the number and location of poles of the prop¬ 
agator and/or intersection points of the running mass. 
Therefore, we decided to repeat the study of TR and BP 
and investigate the propagator in the complex plane in¬ 
cluding all Riemann sheets nearest to the first (physical) 
sheet in order to clarify this problem. The self-energy on 
the unphysical sheet (s) is obtained by analytic continua¬ 
tion. To this end, one first computes the discontinuity of 
the self-energy across the real s-axis, 

DiscII(s)=2* lim ImllAs + ie) , s£R. (7) 
£->o+ t-— 1 

i 

Then, the appropriately continued self-energy n c (s) on 
the next Riemann sheet is obtained via 

n c ( s ) = n( s ) + Discn( s ). (8) 

This expression is valid on the whole Riemann sheet, i.e., 
s is complex-valued. Note that in our case there are three 
thresholds, in successive order corresponding to the de¬ 
cays of ao into irr], KK, and 7 xrj'. These channels will 
be numbered i = 1, 2,3 in the following. Thus, crossing 
the real s-axis at values of s in the interval (sth,i, st/ 1 , 2 ], 
we move from the first to the second Riemann sheet, in 
the following denoted by roman numeral II. Analogously, 
crossing the real s-axis in the interval (sth, 2 , Sth., 3 ], we 
move from the first to the third (III) sheet. Finally, cross¬ 
ing the real s-axis at s > Sth, 3 , we move from the first 
to the sixth (VI) sheet (in the standard notation). Since 
we will also show plots of the spectral function d(x ), we 
recall its definition, 

2x 

d(x) = —— lim ImA(i 2 + ie) , (9) 

7T £->0+ 

where x = y/s. 

B. Spectral functions and poles 

We introduce a dimensionless parameter A £ [0,1] and 
replace the coupling constants in Eq. § by gf -» A gf. In 
consequence, for A = 0 the self-energy vanishes and we 
just obtain the spectral function of the non-interacting 
seed state, i.e., a delta function. The corresponding pole 
lies on the real y^s-axis. Increasing A from zero to 1, the 
interaction is successively increased and we can monitor 
in a controlled manner how the spectral function changes. 


In the following figures, we will show the spectral function 
for the physical value A = 1.0 and for the intermediate 
value A = 0.4. Changing A from zero to one, we will 
also see how the pole of the seed state moves off the real 
axis and other poles emerge, which correspond to the 
dynamically generated resonances. A continuous change 
of A will trace out pole trajectories in the complex y/s- 
plane. The final and physical locations of the poles are 
reached when A = 1.0, which we indicate by a dot in the 
following figures. We consider the three Riemann sheets 
nearest to the physical region (i.e., the first sheet) in one 
figure (a list of the poles corresponding to the resonances 
of interest can be found in Sec. IV). For TR, we use the 
values gi = 1.2952 GeV, g 2 = 0.8094 GeV, g 3 = 0.9461 
GeV, and k 0 = 0.56 GeV, and for BP 5l = 1.7271 GeV, 
g 2 = 1.0975 GeV, g 3 = 1.4478 GeV, and k 0 = 0.56 GeV. 

The results are shown in Fig. [T| We first discuss the 
results for the TR parameterization, and then those for 
BP: 


1. The two panels in the upper row of Fig. [l] show the 
results of TR. For A = 1.0 the spectral function ex¬ 
hibits a narrow peak in the region around 1 GeV 
that was interpreted by TR as the ao(980) reso¬ 
nance. We furthermore observe a broad structure 
above 1.5 GeV. For decreasing coupling strength 
the narrow peak around 1 GeV vanishes, while the 
broad structure becomes more pronounced. It is lo¬ 
cated around 1.4 GeV, which is the location of the 
seed state. The width of the peak decreases with 
A, such that we obtain a delta function for A = 0, 
as expected (not displayed here). 

The behavior described above can also be under¬ 
stood considering the pole structure in the com¬ 
plex y/s-plane. The narrow peak around 1 GeV for 
A = 1.0 corresponds to a pole at s w (1.084 2 — 
i 1.084 • 0.270) GeV 2 , which TR has found on the 
second sheet. This pole is indeed present only if 
A exceeds the critical value A™ « 0.75. The pole 
emerges close to (but not on) the real axis for A™ 
and descends down into the complex plane on the 
second sheet for increasing coupling strength. One 
can interpret this appearance and motion of a pole 
as a feature typical for the kind of dynamical gen¬ 
eration we are interested in. 

However, we also find another pole on the second 
sheet emerging at a large imaginary value of y/s and 
moving up towards the real axis. It first appears 
for A™ ~ 0.84. Its effect on the spectral function 
is hard to discern, since (the absolute value of) its 
imaginary part (ie., its decay width) is still too 
large at A = 1.0. This pole was not reported in 
TR, yet, in Ref. |36j . a similar situation was de¬ 
scribed where the ao(980) was taken to have such a 
behavior, i.e., its pole was coming from the region 
of large negative imaginary parts of y/s and head¬ 
ing towards the real axis. This is, however, not the 
case for the pole of TR, which is dynamically gen- 
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FIG. 1. Spectral functions (left panels) and positions of poles in the complex -^/i-plane (right panels) for the parameter sets 
of TR (upper row) and BP (lower row). Spectral functions are shown for A = 0.4 (dashed grey lines) and A = 1.0 (solid blue 
lines). The pole trajectories of the seed state are indicated by grey dotted or red dashed lines (for details, see the text) and 
the one for the dynamically generated resonance by solid blue lines. The roman numerals indicate the Riemann sheets where 
the respective poles can be found. Final pole positions (A = 1.0) are indicated by solid black dots, pole positions at A Ct i, i.e., 
where the pole i first emerges, are indicated by X. 


erated near the real axis and then shifts towards 
larger (negative) imaginary values of yfs. 

On the third sheet, TR reports another pole. One 
could think that this pole corresponds to the seed 
state, since for A = 0 the pole trajectory starts on 
the real axis at the mass of the seed state. How¬ 
ever, the pole lies on the third sheet, i.e., prior to 
crossing the Trry'-threshold, but its location is above 
that threshold. Therefore, it should not be con¬ 
sidered to induce the broad bump in the spectral 
function. However, there is also a pole on the sixth 
sheet which also starts at the mass of the seed state. 
From its position this pole can also be considered 
to generate the broad resonance shape in the spec¬ 
trum above 1.5 GeV. It is interesting that the pole 
on the third sheet was suggested in TR to corre¬ 
spond to the do (1450) resonance. From our point 
of view, because of the above arguments it is more 


natural to take the pole on the sixth sheet. A close 
inspection of the peak position of the broad bump 
in the spectral function reveals that it corresponds 
more closely to the real part of the pole on the sixth 
sheet than that on the third sheet, which corrobo¬ 
rates our interpretation. 

2. In the lower row of Fig. |T] we present the results 
for the parameter choice in BP. We find that the 
qualitative behavior is very similar to the one in 
TR. Quantitatively, we find that the bump in the 
spectral function corresponding to the ao(980) reso¬ 
nance is now somewhat wider. The broad structure 
at large y/s is more pronounced and now lies around 
2 GeV. For decreasing coupling strength the peak 
becomes narrower and moves towards 1.6 GeV (be¬ 
cause the seed state is located there). 

We find again two poles on the second sheet. The 
right pole appears first for A® ^ ss 0.69 and the left 























5 


one for A®? « 0.66. The parameter set of BP does 
not yield a pole structure from which one can infer 
which pole corresponds to the a 0 (980). Both poles 
give too large widths, and the left one is too light, 
while the right one is too heavy. It seems that both 
of them are relevant in the generation of the bump 
at 1 GeV in the spectrum. Moreover, it does not 
seem to be appropriate to assign the poles on the 
other two sheets to ao(1450). At least within this 
model and with the chosen parameters, the pole 
masses are definitely too high Ea¬ 


rn. SIMPLE EFFECTIVE MODEL WITH 
DERIVATIVE INTERACTIONS 

In the previous section we have re-examined the ap¬ 
proach of TR and BP to dynamically generate resonances 
in the scalar-isovector sector. We now apply the above 
mechanism of dynamical generation of resonances using 
a formulation based on an interaction Lagrangian. 


A. Interactions with derivatives: a lesson from the 
eLSM 

The way a scalar field couples to pseudoscalar states 
depends on the effective approach used. Let us, for in¬ 
stance, consider the coupling of ao to kaons. In chiral 
perturbation theory (chPT) [38], which is based on the 
nonlinear realization of chiral symmetry, only derivative 


couplings of the type K°d^K° can appear in the chi¬ 
ral limit [39]. Away from the chiral limit, a non-derivative 
coupling clqK 0 K 0 appears, too, but its strength is pro¬ 
portional to rn ‘ 2 K , i.e., via the Gell-Mann-Oakes-Renner 
relation proportional to the explicit breaking of chiral 
symmetry by nonzero quark masses. On the other hand, 
if the standard linear sigma model (without vector de¬ 
grees of freedom) is considered, the coupling is only of the 
non-derivative type a[]A" 0 A' 0 . At tree-level both chPT 
and the sigma model can coincide, but when loops are in¬ 
cluded differences arise due to the different s-dependence 
in the amplitudes. 

Studying the spectral function of cf> —► ao (980)7 ~► 
7r°T77 measured by the KLOE Collaboration [40], it was 
shown in Ref. m that a derivative coupling of the type 
a^d^K^d^K 0 seems to be necessary. As we shall demon¬ 
strate below, we come to the same conclusion: a deriva¬ 
tive coupling is necessary for the simultaneous descrip¬ 
tion of both resonances ao(980) and ao(1450). 

Interestingly, an improved version of the linear sigma 
model, called extended Linear Sigma Model (eLSM), nat¬ 
urally contains both non-derivative and derivative cou¬ 
pling terms. This feature is due to the inclusion of (axial- 
hector degrees of freedom in the model, for more details 
see Refs. [20ll2Tj . This model is able to provide a surpris¬ 
ingly good description of the tree-level masses and decay 
widths of meson resonances below 1.7 GeV [120H22j . Fur¬ 
thermore, in this approach the resonance ao(1450) turns 
out to be predominantly a quark-antiquark state with a 
(bare) mass of m ao = 1.363 GeV. 

The Lagrangian for the scalar-isovector sector emerg¬ 
ing from the eLSM has the following form: 


/’eLSM 

U 'a 0 r)TT 

/’eLSM 

^'aQTj'TT 

/-eLSM 

a 0 KK 


= Af SM a°777r° + 3 M a^d„ V d^n 0 + Cf SM 8 ^( 71 °+ rfi*n°) , 

= Af SM a° 0 r]'n 0 + B§ LSM ag0 <1 »/0' 1 tt° + Q LSM ^a[]( n°8^r]' + rfd*n°) , 
= Af SM a%(K°K 0 - K~K+) + Bf SM a^{d II K 0 d l, 'K [) - 8^K~d ll K + ) 
+ CZ LSM d tl a%(K°d ,l K 0 + K°d^K° - K~d fl K + - K + 8^K~) , 


( 10 ) 


where A® LSM , Bf LSM , and (7® LSM are coupling constants 
that are functions of the parameters of the model m- 
Note that both non-derivative and derivative interactions 
appear. The derivatives in front of the fields produce 
an s-dependence in the decay amplitudes, — iA4° LSM (s), 
which enter the tree-level expressions of the decay widths, 


ieLSM/„\ _ ki(s) .. , e LSM^„\|2 


\s) = 


8ns 


-iM^ bM (s)j 2 0(s-s th , t ) , (11) 


which have to be evaluated for s = m 2 o . The amplitudes 
read 

Alf SM (s) = A 


1 1 nc 

~2 Bi 


(s-m^-m^+Cf—s, 
( 12 ) 

where the masses mu , m *2 are the pseudoscalar masses 


in the relevant channels. The parameters of the eLSM 
were determined from a \ 2 -fit to tree-level masses and 
decay widths. So far, no loop corrections were consid¬ 
ered. For a consistent loop calculation one would have to 
perform a new fit of the parameters, which is an inter¬ 
esting project for future work. In any case, one should 
not use the values of the parameters of the eLSM deter¬ 
mined in Ref. |21 in the expressions for A® LSM , i?® LSM , 
and (7® lsm . Therefore, we shall treat the latter as free 
parameters in the following. 

A first attempt to incorporate loop corrections in a 
scheme inspired by the eLSM was presented in Refs. 
HU S3] • There, the s-dependence of the amplitudes was 
completely neglected and a regularization function was 
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introduced, 

- iM? SM (s ) -+ -iM i ( S ) = -iAtf SM K 0 )F j (s). 

(13) 

After that, the imaginary part of the self-energy was com¬ 
puted using the optical theorem, 

Imllj(s) = - v / srf ee (s) = ~ \-iM 2 (s)\ 2 0(s-s t h,i) , 

07 Ty/S 

(14) 

and the real part from the dispersion relation ([2]). As 
shown in Ref. [12] the model yields a width of the seed 
state which is too small. Moreover, no additional pole for 
the ao(980) is dynamically generated. Obviously, neglect¬ 
ing the s-dependence of the amplitudes is an oversimpli¬ 
fication. One has to take into account the derivatives in 
some way; at the same time care is needed when deriva¬ 
tive interactions appear in a Lagrangian, for details see 

App. 0 


B. Effective model with both non-derivative and 
derivative interactions 


We now consider an effective model for the isovector 
states containing the same decay channels as the eLSM 
and including also non-derivative and derivative interac¬ 
tions. The Lagrangian is given by the sum of the follow¬ 
ing terms: 

Ca 0 r, 7T = Aio^tt 0 + B^d^rid^TT 0 , (15) 

£a 0 r,’n = A 2 a^r{ ^ + B 2 a°d^rj'd^n 0 , 

Ca 0 KK = A 3 a° 0 (K°K 0 - K-K+) + B 3 a° 0 (d fl K° 8^ K (> 

- d /t K-d ll K+) . 


Formally it can be obtained by r ewriting the terms pro¬ 
portional to C® LSM in Eq. (llOj) by an integration by 
parts in order to get rid of the derivatives of the a o- 
fields [33]. Subsequently, one replaces the emerging sec¬ 
ond derivatives with the help of the Klein-Gordon equa¬ 
tion, Dtt 0 = —rn^7r° (and similarly for the other pseu¬ 
doscalar fields). Then, Eq. ( fl5| ) gives rise to the following 
s-dependent amplitudes: 


Mf(s) 




Fi(s) , 


(16) 


where we have already included a regularization func¬ 
tion Fi(s) as defined in Sec. [II| We again note that the 
constants A.i and Bi will not be computed from the nu¬ 
merically determined parameters of the eLSM, but will 
be determined in order to produce the masses and de¬ 
cay widths of the resonances under study. Note that in 
chPT the parameters Ai are proportional to the masses 
of the pseudo-Goldstone bosons as A\ oc m\ + m^, 
A 2 oc mjj + m^, , A 3 oc 2 mr K and thus vanish in the chi¬ 
ral limit. Thus, also from this consideration, we expect 
that the derivative terms are sizable and crucial for the 
determination of the resonance poles. 


We computed the real and imaginary part of the self¬ 
energy in two ways. In the first, we applied the method 
outlined in Sec. |III A| i.e., we computed the tree-level de¬ 
cay widths and used the optical theorem from Eq. © 
to obtain the imaginary part of the self-energy. We then 
applied the dispersion relation ([2]) to calculate the corre¬ 
sponding real part. In the second approach, we computed 
the one-loop self-energy directly from the Feynman rules. 
From a comparison, we identified the necessity to intro¬ 
duce subtractions in the first approach, for details see 
APP- 0 Note that the one-loop approximation for the 
self-energy is quite reliable, since vertex corrections can 
be shown to have a negligible effect [35] - 

There are eight parameters in our approach: mo, A = 
V2k 0l and six coupling constants A.i,Bi ( i = 1,2,3). We 
vary the numerical values of m 3 and A within reasonable 
intervals m 0 £ (0.8,1.5) GeV and A € (0.4,1.5) GeV and 
each time perform a fit of the six coupling constants to 
six experimental quantities: one pole in the PDG range 
for ao(980) (in our case y/s = (0.969 — i0.045) GeV) and 
one for a 0 (1450) (in our case y/s = (1.450—iO.135) GeV), 
and the central values of the branching ratios of ao(1450) 
[see Eq. ©]■ By this, all six free parameters can be 
fixed. 

It turns out that there is only a narrow range of suit¬ 
able values of the parameters m 3 and A for which the fit 
of the six coupling constants is possible: approximately 
m 0 £ (0.9,1.2) GeV and A G (0.4,0.9) GeV. Here, ‘ap¬ 
proximately’ refers to the fact that, due to the interde¬ 
pendence of the parameters, the window is not rectan¬ 
gular. However, a small change in m 3 and/or A by 50 
MeV near the borders of the quoted interval does not al¬ 
low one to reproduce the data anymore. Thus, although 
we have eight parameters, we are severely constrained in 
their choice in order to describe the 1=1 resonance. As 
we will see below, the present parameters also explain 
why do (980) couples strongly to kaons. The final values 
for the parameters and coupling constants are: 


m 0 = 1.15 GeV , A 

Ai = 2.52 GeV , B x 

A 2 = 9.27 GeV , B 2 

A 3 = -6.56 GeV , B 3 


= 0.6 GeV , 

(17) 

= -8.07 GeV -1 , 

(18) 

= 9.25 GeV -1 , 


= -1.54 GeV" 1 . 



We rescale these coupling constants by a common factor 
y/\ and compute the corresponding spectral functions. 
The result is shown in the left panel of Fig. [2] We also 
compute the pole trajectories in the complex -^/s-plane 
by varying A from zero to one. The following comments 
are in order: 


1. The spectral function shows a narrow peak for 
A = 1.0 at a value of x = y/s slightly smaller than 1 
GeV, which can be interpreted as the ao(980). The 
form is distorted by the nearby A'A'-threshold and 
resembles the Flatte distribution [3E], see also Refs. 
Bung and refs, therein. The pole corresponding 
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0.5 1.0 1.5 2.0 

x [GeV] 



FIG. 2. In the left panel we show the spectral functions for three different values of A. In the right panel we display pole 
trajectories obtained by varying A from zero to one. Black dots indicate the position of the poles for A = 1.0. The X indicates 
the pole position for A c , i.e., when the pole first emerges. The roman numeral indicates on which sheet the respective pole can 
be found. 


to this peak lies on the second sheet and has coor¬ 
dinates 


y/s = (0.970 - 10.045) GeV , (19) 

i.e., we find the ao(980) to have a mass of TO po i e = 
0.969 GeV and a width of r po i e = 0.090 GeV. This 
pole appears only if A exceeds A c ~ 0.52 (note that 
the pole trajectory is very different from the one 
reported in Ref. ED- The corresponding position 
is indicated by an X in the right panel of Fig. [2} The 
important thing here is that, in contrast to what we 
have found for the TR and BP parametrizations, 
there is only one pole for the ao(980), and thus no 
ambiguity which one should be identified with this 
resonance [49] . 


The present study demonstrates that, by starting with 
a unique seed state, it is indeed possible to find two 
poles for the isovector states, both of which reproduce 
the masses and widths of do(980) and ao(1450) reason¬ 
ably well. 


C. Branching ratios and coupling constants for 

do (980) 


For completeness, we report the branching ratios of our 
effective model by using the tree-level decay widths ob¬ 
tained from the optical theorem (14). The partial widths 


are evaluated at the peak value of the spectral function 
above 1 GeV, m pea k = 1.419 GeV. For the resonance 
do (1450) this leads to 


2. There is also a broad structure around 1.5 GeV 
which corresponds to the resonance do (1450). For 
decreasing A, both peaks merge and settle around 
1.15 GeV, where the seed state is located. 

3. As expected, there is (only) one pole present on the 
third sheet with coordinates yfs = (1.400 — i0.141) 
GeV. However, as in TR and BP, we find a pole on 
the sixth sheet, too. Its coordinates are 

7s = (1.456 - *0.134) GeV , (20) 

or TTipoie = 1.450 GeV and r po i e = 0.270 GeV. This 
is the pole which is responsible for the peak around 
1.5 GeV in the spectrum, and thus we assign it to 
the do (1450). However, since the pole on the third 
sheet also reproduces mass and width of do (1450) to 
reasonable accuracy, it is in principle possible to re¬ 
gard this one as the pole corresponding to a o (1450), 
too. 


ptree 

1 a 0 ->7)'7r 

ptree 
1 ao—>rjir 


~ 0.44 , 


ptree 

ap^KK 

ptree 

ao^rjiv 


~ 0.96 , 


( 21 ) 


which can be compared to the experimental values [I]: 
Ta °-+i'” = 0.35 ±0.16 , T “°-* kk = 0.88 ±0.23 . (22) 


- CLQ—>r)7T 


- aQ—tij-n 


Concerning the resonance do (980), we give the follow¬ 
ing estimates for the coupling constants in the nr]- and 
A'R-channels: We calculate the partial widths r* ree (s), 
this time with y/s equal to the peak mass of the spectral 
function below 1 GeV, m pea k = 0.956 GeV. Then, Eq. 
(14) is used to solve for the absolute values of the ampli¬ 


tudes. The result is multiplied with the root of the wave- 
function renormalization factor, \[Z = 0.652, which is 
its value at the Breit-Wigner mass of the do (980). Thus 
we obtain the coupling constants in the nr]- and KK- 
channels as 


3 ^ = 2.496 GeV , g Kit = 6.012 GeV 


(23) 

























It is remarkable that the coupling of do (980) to kaons 
turns out to be sizably larger than the coupling to nr]. 
This is in agreement with various other works on this 
topic HE EM EH]: Virtual kaon-kaon pairs near the kaon- 
kaon threshold are important for the dynamical genera¬ 
tion of the resonance ao(980). 


IV. CONCLUSIONS 

Experimental data exhibit several puzzling facts about 
the light scalar mesons: / 0 (500) (or cr) and Kg (800) have 
large decay widths, while / o (980) and do(980) are nar¬ 
row but their spectral functions show threshold distor¬ 
tions due to the nearby A'A-threshold. It is nowadays 
recognized that these states do not fit into the ordinary 
qq picture based on a simple representation of SU( 3)- 
flavour symmetry [50], yet there is no consensus on the 
precise mechanism which generates them. One can also 
regard these states as four-quark objects, for example as 
tetraquarks pi m EH E3 EH] or as dynamically gener¬ 
ated states. The latter are states which are not present in 
the original formulation of a hadronic model but appear 
when calculating loop corrections mmmm- Indeed, 
the interpretation of light scalar states as loosely bound 
molecular states |29| is also an example of dynamical gen¬ 
eration. [For another interpretation of light scalar states, 
see e.g. Ref. [51].] 

A particular type of dynamical generation is that of 
‘image’ or ‘companion poles’. We have concentrated on 
such a method in this work and have applied it to the 
scalar-isovector sector. Our results demonstrate that it 
is in fact possible to correctly describe the resonances 
do(980) and do(1450) in a unique framework, where orig¬ 
inally only a single quark-antiquark seed state is present. 

Besides that, we have also repeated the previous cal¬ 
culations of Tornqvist and Roos m and Boglione and 
Pennigton m These studies have been extended by us 
to the complex plane on all Riemann sheets nearest to the 
first, physical sheet. A summary of our results and, for 
comparison, those of Tornqvist and Roos, and Boglione 
and Pennington, can be found in Tab. |T] Our results are 
based on an effective Lagrangian approach that includes 
both derivative and non-derivative interaction terms, see 
Eq. |I5| ), inspired by the extended Linear Sigma Model 
(eLSM), and show that both terms are necessary and 
equally important [20, El] ■ 

Note that the formulation of dynamical generation ap¬ 
plied here is related, but not equal to the one described 
in Ref. [5]. In the latter the scattering amplitude is com¬ 
puted from an effective Lagrangian (derived from chi¬ 
ral perturbation theory and containing only pseudoscalar 
mesons) and then unitarized; this process of unitariza- 
tion generates, for instance, the pole of the ao(980) in 
the scalar-isovector sector. Yet, it is a priori not pos¬ 
sible to know if the resulting state is in fact a quark- 
antiquark or a four-quark resonance, and if it can be 
linked to the heavier do(1450) state or not; see Ref. [17: 


for a detailed discussion of this issue. However, further 
studies within this scheme were performed in Ref. i by 
including an octet of bare resonances masses around 1.4 
GeV. It was found that the physical do (1450) in fact orig¬ 
inates from this octet, giving a clear statement about its 
nature, which is in agreement with our results. 

In this work we have concentrated on the isovector 
sector 1=1. However, the very same mechanism is ap¬ 
plicable for the other light and heavy scalar states. As a 
consequence, all light scalars are dynamically generated 
states. In particular, it seems promising to extend the 
present study in the low-energy regime into two direc¬ 
tions: (i) The isodoublet, i.e., by describing the reso¬ 
nances ATq( 800) and A)) (1430) in a similar framework. 
The pole of A~q( 800) is not yet very well known and 
there is need of improved analyses from different direc¬ 
tions. (ii) The scalar-isoscalar sector, where the reso¬ 
nances /o(500) and /o(980) should be dynamically gen¬ 
erated. In this case, /o(1370), /o(1500), and /o(1710) 
would be predominantly a non-strange quarkonium, a 
strange quarkonium, and a scalar glueball, respectively. 

Another interesting subject is the study of dynami¬ 
cal generation in the framework of puzzling resonances 
in the charmonium sector [52] . see for example Ref. [53] 
and refs, therein. Namely, a whole class of mesons, called 
X , Y , and Z states, has been experimentally discovered 
but is so far not fully understood [M, [35]. As demon¬ 
strated in Ref. [SB] for the case of X (3872), some of the X 
and Y states could emerge as companion poles of quark- 
antiquark states. 
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Appendix A 

In this Appendix, we compute the one-loop self-energy 
in the case of derivative interactions. We first derive the 
interacting part of the Hamiltonian from the Lagrangian 
via a Legendre transformation. We shall see that the 
derivative interactions give rise to new interaction ver¬ 
tices. We demonstrate that, in a perturbative calculation 
of the one-loop self-energy, these terms are necessary to 
cancel additional terms arising from contractions of gra¬ 
dients of fields. This proves that, at least at the one-loop 
level, it is justified to apply standard Feynman rules with 
the derivative interaction in £i nt . We shall also demon¬ 
strate that a computation of the self-energy via the dis¬ 
persion relation § may require subtraction constants to 
agree with the perturbative calculation using Feynman 
rules. 
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m po le [GeV] 

ao(980) 

Tpoie [GeV] 

a 0 (1450) 

m P oie [GeV] 

r po ie [GeV] 

TR HU 

1.084 

0.270 

1.566 

0.578 

bp U3 

1.186* 

0.373* 

1.896 

0.250 

Our results 

0.969 

0.090 

1.450 

0.270 

PDG [J 

0.980 ±0.020 

0.050 to 0.100 

1.474 ±0.019 

0.265 ±0.013 

*In order to compare to TR, the right pole 

on the second sheet was chosen. 




TABLE I. Numerical results for the pole coordinates in the scalar-isovector sector in TR, BP, and our effective model, compared 
to the PDG values. In the case of the ao(1450), the poles listed for TR and BP are located on the third sheet, while our pole 
lies on the sixth sheet. All poles for the ao(980) are found on the second sheet. Note that all poles listed for BP were obtained 
performing the analytic continuation of the propagator given by BP. 


Canonical quantization 

Let us consider a theory with two scalar fields, S and 
4>, which allows for the decay process S —> $<j). Conse¬ 
quently, the Lagrangian is 

£ = £s + + Ant > (Ai) 

where 

C S = \ (duSWS - M 2 S 2 ) , 

£<t> = \ {d l i<t>d li <t> - m 2 cj) 2 ) , 

Ant = gSd^d^ . (A2) 

For perturbative calculations of S'-matrix elements or 
Green’s functions, however, one needs the interaction 
part of the Hamilton operator in the interaction picture. 
We derive this operator via a Legendre transformation of 
C and subsequent canonical quantization in the interac¬ 
tion picture. As a byproduct of this calculation we will 
explicitly show that the derivative interactions invalidate 
the commonly used relation "Hint = — Ant m- 
The canonically conjugate fields are 

* s = em =a ° s - (A3) 

f\r 

= d 0 <j> + 2gSd°(j) = (1 + 2gS)d°(j) . 

The Hamiltonian is defined via a Legendre transforma¬ 
tion of C, 

H = 7r sd°S + 7r^<9°</> — C (A4) 

= Kstts + * VS ■ VS + \m 2 S 2 + ^7T 0 7T 0 (1 + 2 gS)- 1 

+ ^Vtp • Vcj) + ^ m 2 (f i 2 + gSV<j> ■ V(j> . 

For a perturbative calculation, we need to expand the 
denominator (1 ±2<?iS')~ 1 and obtain the interaction part 
of the Hamiltonian as 

■Hint = -gSn,p 7T0 + gSV(j) ■ V(j) + 2g 2 S 2 n+ 0(g 3 ) . 

(A5) 


We may now quantize in the Heisenberg picture (indi¬ 
cated by a superscript H at the respective operators). 
This is commonly done by promoting fields to operators 
S -A S H , <f> —> (j) H , tts —> , 7T0 —> ir^, and postulating 
certain commutation relations for these operators. How¬ 
ever, in perturbation theory we need the operators in the 
interaction picture. The following relations hold, 

S 1 = US H Ui , ft = Uj) H U^ , (A6) 

*S = forff/t , tt' = Utftf , 

where U = e lHot e - lHt , g ^g time-evolution operator 
that relates operators in the Heisenberg picture with 
those in the interaction picture. Finally, replacing = 
d 0 ^ 1 , 7= d°ft this results in 

Kt = -£fnt + 2 g 2 S I S I d 0 ftd°ft + <D(g 3 ) . (A7) 

As advertised, the second term spoils the standard re¬ 
lation 7-Li n t = —Tint • This term corresponds to a four- 
point vertex, so it will not appear in the tree-level decay 
S —► 4>4>. In contrast, in the one-loop self-energy, it will 
give rise to an additional tadpole contribution. 


Perturbative calculation of the one-loop self-energy 


We now turn to the self-energy n(s) of the field S. At 
one-loop level, the Feynman rules applied to Hint tell us 
that we will have two contributions. The first contribu¬ 
tion comes from taking two three-point vertices of Ant 
where the <f> legs are joined in a manner which gives a 
1PI diagram. A covariant derivative acts on each (f> leg 
at each vertex. The second contribution is a tadpole term 


arising from the four-point vertex in Eq. (A7), which has 


two time derivatives on the internal leg. This can be 
graphically depicted as follows: 




— ~(£**) + 


(A8) 
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The usual Feynman propagator is defined as a contrac¬ 
tion of two fields: 


( j>{xi)(j) (x2) = (O|T{0(xi)<£(x 2 )}|O) 

= (O|0(xi)0(x 2 )|O)0(xi - x°) 

+ (O|0(x 2 )^(xi)|O)0(x° - x°) 
= iAp( x i “ x 2) 


= i 


d 4 p e -ip-(x 1 -x 2 ) 

(27r) 4 p 2 — m 2 + ie 


(A9) 


However, in the tadpole diagram, we have the contraction 
of two fields on each of which acts a time derivative. 
Because time-ordering has no effect at the same space- 
time point, we obtain: 

(0\T{d$4>(x)d°’ x <j>{x)}\0) = (0\dS$(x)dP'*Ux)\0) 


= l 


j4 

d p 


E l 


(27t) 4 p 2 — m 2 + ie 

(A10) 


In order to obtain this result we inserted the standard 
Fourier decomposition of the field operators 


(j)(x) = 


d 3 p 1 


(«P e 


+ a D e 


c ), (All) 


where p° = E p = ^/p 2 + m 2 is the on-shell energy. Thus, 
a time derivative acting on a field operator brings down a 
factor of ±i times the on-shell energy in the correspond¬ 
ing Fourier representation. 

The result (A10) is, however, identical if we just act 


with the time derivatives on the standard Feynman prop¬ 
agator (A9): 


5o9°’ x (O|0(x)0(x)|O) = lim d^d°’ X2 (0|T{^(xi)<(>(x 2 )}|0) 

Xl—>X2 


= l 


i4 

d p 


E l 


(27t) 4 p 2 — to 2 + ie 


(A12) 


In order to prove this, it is convenient to first perform 
the po integration in Eq. (A9) and then take the time 
derivatives. The equivalence of Eqs. (A10) and (A12) is 
graphically depicted as 


<0|T{«x)d°’^(x)}|0> 



(A13) 



In the perturbative series of the full propagator of the 
S-field, this tadpole contribution appears in combination 
with two free S'-field propagators (where we omit the su¬ 
perscript S ): 


—■ 2 J da;' iA(xi — x')(0|T{<9o 4 >(x')d°’ x 0(x , )}|O}iA(x / — x 2 ) (A14) 

= i2g 2 ■ 2 J da/ A(xi — x')A(x' — x 2 )(0|T{c?o <t>(x')d 0,x (/>(x')}|0) . 


The factor —2 ig 2 is the factor accompanying the four-point vertex, cf. Eq. (A7). A factor of two arises because each 
S propagator can be joined with either one of the S legs at the vertex. 


We now compute the first diagram in Eq. (A 81 . To this end, we need contractions of gradients of the ^-fields. These 


can be expressed in terms of gradients acting on the standard Feynman propagator. The gradient of the Feynman 
propagator (A9| is 


i> 9 X2 A^(xi - x 2 ) = 5 X2 (O|T{^(xi)0(x 2 )}|O) 

= (0|T{</>(xi)<9 X2 ^(x 2 )}|0) - r] v0 5{x\ - x§)(0| ),</>(x 2 )] |0) , 


(A15) 


=o 


where we used the explicit definition of the time-ordered product. The last term vanishes on account of the delta 
function, since it is an equal-time commutator of two <(>-fields m- Taking another gradient leads to 

id^ 1 d* 2 Ap(xi — x 2 ) = 9 Xl (O|T{<(>(x 1 )9 X2 0(x 2 )}|O) (A16) 

= (0|T{<9 Xl <(>(xi)<9 X2 ^(x 2 )}|0) + ri^, 0 6(xi - x°)(0| [0(xi), 5 X2 0(x 2 )] |0) . 

"-V-' 

7^0 


Now the second term does not vanish if v = 0, because then it involves a commutator of a field with its canonically 
conjugate field El. 


VpO^( x i — ®§)(0| [<^(xi), 9 X2 <()(x 2 )] |0) = irivoVvoSW{x! - x 2 ) . 


(A17) 
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Collecting terms, we can express the contraction of two gradients of the 0-field as 

(0|7~0(^i)<9^ 2 0 (^ 2 )}|0) = id Xl d X2 A^(xi - x 2 ) - *?7 M oWo6 (4) (xi - x 2 ) ■ 


(A18) 


In the perturbative series of the full propagator for the S'-field, the first diagram in Eq. (A8) also appears in the 
combination with two 0-field propagators: 


tig? 

2 ! 


■2-2 / dx' / dx" iA(x 1 -x l * * * * * * * )iA{x"-x 2 )(0\T[d*(l)(x , )df(l)(x")}\0){0\T{d IJ ’’ x '(l)(x')d , '’ x "<l){x")}\0). (A19) 


Two factors of — ig originate from the three-point vertices in £i nt . The factor of 1/2! arises because the diagram is 
second order in perturbation theory. A factor of two arises because each S propagator can be joined with either one 
of the S legs at the vertex. Finally, another factor of two comes from the fact that the two 0 lines at one vertex 
can be joined with corresponding lines at the other vertex in two different ways. Successively inserting Eq. (A18) we 
compute 

20 J da/ J da/ 7 A(xi — x 7 )A(x" — x 2 ) id x d x Ap(x' — x") — ir) tl oVv 0 ^ 4 '\x' — x") (0|T{$ /7,x 0(x 7 )d ! '’ x 0(x 77 )}|O) 

= 20 J dx' J dx" A(aq - a/) A (a/ 7 - x 2 ) id x> df A*(x' - x 77 )(0|T{d /J ’ x >(x 7 )d I/ ’ x >(x 77 )}|0) 

— 120 J dx' A(aq — x 7 )A(x 7 — x 2 )(0| T{d 0,x (f)(x')d°’ x 0(x 7 )}|O) 

= -20 J dx 7 J dx" A(aq - x 7 )A(x" - x 2 ) d£df A*(x' - x”)d^ x 'd v ’ x "A*(x' - x") 

— i2g 2 J d x ' J dx 77 A(xi — x 7 )A(x 77 — x 2 )idy d$ Ap(x' — x 77 )0 4 )(x 7 — x 77 ) 

— i2g 2 J dx 7 A(xi — x 7 )A(x 7 — x 2 ){Q\T{d° lX 0(x 7 )0°’ x 0(x 7 )}|O) . (A20) 


With Eq. (A13) one realizes that the last two terms are identical. The final result is 

—20 J dx' J dx" A(xi - x')A{x" - x 2 )d x 'd x " Ap(x' - x")d^ x 'd v ' x " A%{x' - x") 
— i2g 2 • 2 J dx 7 A(xi — x 7 )A(x 7 — x 2 )(0|T{i9 0 ’ x 0(x 7 )<9 o,x 0(x 7 )} 10) , 


(A21) 


which can be graphically depicted as 


“(S) @- 


d.dj 




/d^d" 


d 0 a° 


(A22) 

Obviously, the second diagram cancels the tadpole con¬ 
tribution, Eq. (A13), in the one-loop self-energy from Eq. 

^ 81 - 


In summary, a derivative interaction in £i nt produces 
an additional term in the interaction Hamiltonian and 
thus, after quantization, an additional vertex which has 
to be taken into account in perturbative calculations via 
Feynman rules. In the one-loop self-energy, this vertex 
leads to a tadpole diagram. Nevertheless, carefully com¬ 
puting contractions between gradients of the field opera¬ 
tors we demonstrated that these lead to a term which ex¬ 
actly cancels the tadpole diagram. The remaining contri¬ 
bution is exactly equal to the self-energy when computed 


I 

with standard Feynman rules using £; nt and derivatives 

acting on the usual Feynman propagators. 

We did not deliver a rigorous proof of this tadpole can¬ 

cellation to all orders in perturbation theory. However, 

since this seems to be just a demonstration of the validity 
of Matthews’s theorem [55] which was investigated e.g. 
in Refs. [55H62| . we also expect a similar cancellation to 

work at higher orders in perturbation theory. 


One-loop self-energy from a dispersion relation 

The second way to compute the self-energy is via the 
dispersion relation |2]). To this end, one needs the imag¬ 
inary part of the self-energy in order to compute the real 
part. The imaginary part can be inferred from the de¬ 
cay width through the optical theorem. For the one-loop 
self-energy, the cutting rules imply that the decay width 
needs only to be known at tree-level: 
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/ 


dr 



— Im 






(A23) 


The second equality arises from Eq. ( |A22| ) and the fact that the tadpole has no imaginary part. 

The calculation of the tree-level decay width in momentum space proceeds by replacing derivatives d L —> (the 

lower/upper sign stands for incoming/outgoing particles) in the Lagrangian (A2), z.e., in our simple model the decay 
amplitude reads 



= ~Pl'P2 


(A24) 


The factor two arises from the two identical particles in 
the outgoing channel. The blob in the left diagram repre¬ 
sents the vertex as given by Eq. (A71, while in the middle 
diagram the above replacement was performed in order 
to calculate the expression on the right-hand side. Note 
that the factor (s — 2m 1 2 3 4 5 6 7 )/2 appears in a similar form in 
Eq. (12) with coupling constant g = 5? LSM . Neverthe¬ 
less, since s = M 2 , the tree-level decay width is simply a 

constant. _ 

Returning now to the imaginary part (A23) of the self¬ 
energy, we observe that, on account of the fact that the 
tadpole does not contribute to the imaginary part, with 
the dispersion relation® one actually only computes the 
first diagram in Eq. (A8), but misses the tadpole contri¬ 
bution. In other words, as we have demonstrated above, 
the first diagram in Eq. (A8) contains precisely the tad¬ 
pole contribution, but with opposite sign, cf. Eq. (A22). 
Consequently, we need to add this tadpole to the (real 
part of the) self-energy as computed via the dispersion 


relation, in order to have the latter agree with the result 
obtained from the perturbative calculation. We remark 
in passing that the emergence of a tadpole can also be 
explicitly demonstrated by cleverly manipulating the ex¬ 
pression for II(s) as it results from the Feynman rules. 
In the case of our effective model from Sec. IIIII we com¬ 
puted the self-energies precisely in the manner explained 
above, i.e., from a dispersion relation and adding the 
corresponding tadpole diagrams. 


We conclude by remarking that, in the eLSM La¬ 
grangian for the scalar-isovector state, derivatives also 
occur in front of the decaying field ao; thes e ar e the 
terms with coupling constants Cf LSM in Eq. (10). All 
that has been stated above applies also in this case, with 
the exception that, besides constant tadpole terms, also 
s-dependent contributions appear in the expression for 
II(s). We note that all these additional contributions 
also cancel in a similar way as we discussed above. 


[1] K. A. Olive et al. (Particle Data Group), Chin. Phys. C 
38, 090001 (2014). 

[2] C. Amsler and N. A. Tornqvist, Phys. Rep. 389, 61 
(2004); E. Klempt and A. Zaitsev, ibid. 454, 1 (2007), 
arXiv:0708.4016 [hep-ph] 

[3] F. E. Close and N. A. Tornqvist, J. Phys. G 28, R249 
(2002), arXiv:hep-ph/0204205 

[4] J. R. Pelaez, (2015), arXiv: 1510.00653 [hep-ph] J. R. 
Pelaez and G. Rios, Phys. Rev. Lett. 97, 242002 (2006), 
arXiv:hep-ph/0610397; J. R. Pelaez, ibid. 92, 102001 
(2004), arXiv:hep-ph/0309292 

[5] J. A. Oiler and E. Oset, Nucl. Phys. A 620, 438 (1997), 
arXiv:hep-ph/9702314; J. A. Oiler, E. Oset, and J. R. 
Pelaez, Phys. Rev. Lett. 80, 3452 (1998), arXiv:hep- 
ph/9803242 Phys. Rev. D 59, 074001 (1999), [Erratum- 
ibid. D60, 099906 (1999); Erratum-ibid. D75, 099903 
(2007)], arXiv:hep-ph/9804209 

[6] J. A. Oiler and E. Oset, Phys. Rev. D 60, 074023 (1999), 
arXiv:hep-ph /9809337 

[7] M. Jamin, J. A. Oiler, and A. Pich, Nucl. Phys. B 
587, 331 (2000), arXiv:hep-ph/0006045; M. Albaladejo 
and J. A. Oiler, Phys. Rev. Lett. 101, 252002 (2008), 


arXiv:0801.4929 [hep-ph] 

[8] E. van Beveren, D. V. Bugg, F. Kleefeld, and G. Rupp, 
Phys. Lett. B 641, 265 (2006), arXiv:hep-ph/0606022 

[9] D. Morgan and M. R. Pennington, Phys. Rev. D 48, 1185 
(1993). 

[10] E. van Beveren, T. A. Rijken, K. Metzger, C. Dullemond, 
G. Rupp, and J. E. Ribeiro, Z. Phys. C 30, 615 (1986), 
arXiv:0710.4067 [hep-ph] 

[11] N. A. Tornqvist, Z. Phys. C 68, 647 (1995), arXiv:hep- 
ph/9504372; N. A. Tornqvist and M. Roos, Phys. Rev. 
Lett. 76, 1575 (1996), arXiv:hep-ph/9511210vl| 

[12] M. Boglione and M. R. Pennington, Phys. Rev. Lett. 79, 
1998 (1997), arXiv:hep-ph/9703257 

[13] M. Boglione and M. R. Pennington, Phys. Rev. D 65, 
114010 (2002), arXiv:hep-ph/0203149 

[14] Z.-Y. Zhou and Z. Xiao, Phys. Rev. D 83, 014010 (2011), 
arXiv: 1007.2072 [hep-ph] 

[15] Z.-H. Guo and J. A. Oiler, (2015), arXiv: 1508.06400 
[hep-ph] 

[16] M. Albaladejo and J. A. Oiler, Phys. Rev. D 86, 034003 
(2012), arXiv:1205.6606 [hep-ph] 




























13 


[17] F. Giacosa, Phys. Rev. D 80, 074028 (2009), 

arXiv:0903.4481 [hep-ph] 

[18] C. Amsler and F. E. Close, Phys. Lett. B 353, 385 (1995); 
W. Lee and D. Weingarten, Phys. Rev. D 61, 014015 
(1999), arXiv:hep-lat/9910008; F. E. Close and A. Kirk, 
Eur. Phys. J. B 21, 531 (2001), arXiv:hep-ph/0103173 
H.-Y. Cheng, C.-K. Chua, and K.-F. Liu, Phys. Rev. 
D 74, 094005 (2006), arXiv:hep-ph/0607206i F. Gia¬ 
cosa, T. Gutsche, V. E. Lyubovitskij, and A. Faessler, 
72, 094006 (2005), arXiv:hep-ph/0509247i F. Giacosa, 
T. Gutsche, and A. Faessler, Phys. Rev. C 71, 025202 
(2005), arXiv:hep-ph/0408085 

[19] F. Giacosa, Phys. Rev. D 75, 054007 (2007), arXiv:hep- 
ph/0611388 

[20] S. Gallas, F. Giacosa, and D. H. Rischke, Phys. Rev. D 
82, 014004 (2010), arXiv:0907.5084 [liep-ph]; D. Pargan- 
lija, F. Giacosa, and D. H. Rischke, 82, 054024 (2010), 
arXiv:1003.4934 [hep-ph]; S. Janowski, D. Parganlija, 
F. Giacosa, and D. H. Rischke, 84, 054007 (2011), 
arXiv:1103.3238 [hep-ph] 

[21] D. Parganlija, P. Kovacs, G. Wolf, F. Giacosa, and 
D. H. Rischke, Phys. Rev. D 87, 014011 (2013), 
arXiv:1208.0585 [hep-ph] 

[22] S. Janowski, F. Giacosa, and D. H. Rischke, Phys. Rev. 
D 90, 114005 (2014), arXiv: 1408.4921 [hep-ph] 

[23] L.-C. Gui, Y. Chen, G. Li, C. Liu, Y.-B. Liu, J.-P. Ma, 
Y.-B. Yang, and J.-B. Zhang, Phys. Rev. Lett. 110, 
021601 (2013), arXiv: 1206.0125 [hep-lat] 

[24] R. L. Jaffe, Phys. Rev. D 15, 267 (1977); 15, 281 (1977). 

[25] L. Maiani, F. Piccinini, A. D. Polosa, and V. Ri- 
quer, Phys. Rev. Lett. 93, 212002 (2004), arXiv:hep- 
ph/0407017 

[26] G. ’t Hooft, G. Isidori, L. Maiani, A. D. Polosa, 
and V. Riquer, Phys. Lett. B 662, 424 (2008), 
arXiv:0801.2288 [hep-ph] 

[27] A. H. Fariborz, Int. J. Mod. Phys. A 19, 2095 (2004), 
arXiv:hep-ph/0302133, A. H. Fariborz, R. Jora, and 
J. Schechter, Phys. Rev. D 72, 034001 (2005), arXiv:hep- 
ph/0506170 

[28] M. Napsuciale and S. Rodriguez, Phys. Rev. D 70, 
094043 (2004), arXiv:hep-ph/0407037 

[29] O. Krehl and J. Speth, Nucl. Phys. A 623, 162c (1997); 
C. M. B. Voloshin and L. B. Okun, Pisma Zh. Eksp. 
Teor. Fiz. 23, 369 (1976), [JETP Lett. 23, 333 (1976)]; 
V. Baru, J. Haidenbauer, C. Hanhart, Y. Kalashnikova, 
and A. E. Kudryavtsev, Phys. Lett. B 586, 53 (2004), 
arXiv:hep-ph/0308129, T. Branz, T. Gutsche, and 
V. E. Lyubovitskij, Phys. Rev. D 78, 114004 (2008), 
arXiv:0808.0705 [hep-ph]; S. Krewald, R. H. Lem- 
mer, and F. P. Sassen, 69, 016003 (2004), arXiv:hep- 
ph/0307288 

[30] R. E. Peierls, in Proceedings of the Glasgow Conference 
on Nuclear and Meson Physics (Pergamon Press London, 
1954). 

[31] N. A. Tornqvist and A. D. Polosa, Nucl. Phys. A 692, 
259c (2001), arXiv:hep-ph/0011109 

[32] Note that we use a different sign convention for the prop¬ 
agator A(s) and the self-energy II(s) than Refs. [Ill |T3j . 
In the first reference the scattering amplitude is studied, 
but only its denominator is important here, which equals 
the inverse propagator. 

[33] M. Harada, F. Sannino, and J. Schechter, Phys. Rev. 
Lett. 78, 1603 (1997), arXiv:hep-ph/9609428; N. A. 
Tornqvist and M. Roos, 78, 1604 (1997), arXiv:hep- 


ph/9610527 

[34] TR and BP did not quote values for mu, m; 2 • In this 
work, we consistently used the isospin-averaged numeri¬ 
cal values given in the PDG from 2002, the year of publi¬ 
cation of BP. Note that these values differ from the ones 
used in TR, so that our results for the pole positions 
slightly differ numerically from theirs. 

[35] We apply the usual parameterization for propagator 

poles, Spole — ^'pole i TTlpole 1 pole 

[36] E. van Beveren and G. Rupp, Eur. Phys. J. C 10, 469 
(1999), arXiv:hep-ph/9806246 

[37] From the discussion of the poles, we can furthermore con¬ 
clude that crossings of the running mass are not really 
indicative of poles in the propagator. BP reports three 
crossings, the first at a mass value close to 1 GeV, the 
second one around 1.4 GeV, while a third one is located 
around 1.8 GeV. The latter was discarded in BP as un¬ 
physical [see Ref. El for more details]. From our point 
of view it is not possible to unambiguously assign poles 
to these crossings. 

[38] J. Gasser and H. Leutwyler, Ann. Phys. 158, 142 (1984); 
see also S. Scherer, Adv. Nucl. Phys. 27, 277 (2003), 
arXiv:hep-ph/0210398 and refs, therein 

[39] G. Ecker, J. Gasser, A. Pich, and E. D. Rafael, Nucl. 
Phys. B 321, 311 (1989). 

[40] F. Ambrosino et al. (KLOE Collaboration), (2007), 
arXiv:0707.4609 [hep-ex] 

[41] F. Giacosa and G. Pagliara, Nucl. Phys. A 812, 125 
(2008), arXiv:0804.1572 [hep-ph]; 833, 138 (2010), 
arXiv:0905.3706 [hep-ph] 

[42] T. Wolkanowski, Acta Phys. Polon. B Proceed. Suppl. 8, 
273 (2015), arXiv:1410.7022 [hep-ph], 

[43] T. Wolkanowski and F. Giacosa, Acta Phys. Polon. B 
Proceed. Suppl. 7, 469 (2014), arXiv:1404.5758 [hep-ph] 

[44] Quite remarkably, terms with derivatives acting on the 
decaying field are pretty peculiar, see App. [A] for some 
remarks. 

[45] J. Schneitzer, T. Wolkanowski, and F. Giacosa, Nucl. 
Phys. B 888, 287 (2014), arXiv:1407.7414 [hep-ph] 

[46] S. M. Flatte, Phys. Lett. B 63, 224 (1976); 63, 228 

(1976). 

[47] F. Giacosa and G. Pagliara, Phys. Rev. C 76, 065204 
(2007), arXiv:0707.3594 [hep-ph] 

[48] V. Baru, J. Haidenbauer, C. Hanhart, A. E. Kudryavt¬ 
sev, and U.-G. Meifiner, Eur. Phys. J. A 23, 523 (2005), 
arXiv:nucl-th/0410099 

[49] There are two additional poles in the relevant part of the 
complex plane which will not be displayed and discussed 
here: ( i ) a pole deep in the imaginary region on the sec¬ 
ond sheet, and (ii) a pole close to the imaginary axis on 
the sixth sheet. Both have no physical impact. Quite in¬ 
terestingly, we do not observe any virtual bound states. 
Such poles were described by two of us within theories 
without derivative interactions, e.g. in Refs. |63l |64| . 

[50] S. Godfrey and N. Isgur, Phys. Rev. D 32, 189 (1985). 

[51] P. Minkowski and W. Ochs, Eur. Phys. J. C 39, 71 
(2005), arXiv:hep-ph/0404194, Nucl. Phys. Proc. Suppl. 
121, 123 (2003), arXiv:hep-ph/0209225; Eur. Phys. J. C 
9, 283 (1999), arXiv:hep-ph/9811518 

[52] M. R. Pennington and D. J. Wilson, Phys. Rev. D 76, 
077502 (2007), arXiv:0704.3384 [hep-ph] 

[53] N. Brambilla et al, Eur. Phys. J. C 71, 1534 (2011), 
arXiv: 1010.5827 [hep-ph] 







14 


[54] E. Braaten, C. Langmack, and D. H. Smith, Phys. Rev. 
D 90, 014044 (2014), arXiv:1402.0438 [hep-ph]; G. T. 
Bodwin et al, (2013), arXiv:1307.7425 [hep-ph] 

[55] L. Maiani, A. D. Polosa, and V. Riquer, Phys. Rev. Lett. 
99, 182003 (2007), arXiv:0707.3354 [hep-ph] 

[56] S. Coito, G. Rupp, and E. van Beveren, Eur. Phys. J. C 
73, 2351 (2013), arXiv: 1212.0648 [hep-ph] 

[57] W. Greiner and J. Reinhardt, Field Quantization 
(Springer-Verlag Berlin Heidelberg, 1996). 

[58] P. T. Matthews, Phys. Rev. 76, 684 (1949). 


[59] C. Bernard and A. Duncan, Phys. Rev. D 11 , 848 (1975). 

[60] D. Barua and S. N. Gupta, Phys. Rev. D 16, 413 (1977). 

[61] J. Z. Simon, Phys. Rev. D 41, 3720 (1990). 

[62] C. Grosse-Knetter, Phys. Rev. D 49, 6709 (1994), 
arXiv:hep-ph/9306321 

[63] F. Giacosa and T. Wolkanowski, Mod. Phys. Lett. A 27, 
1250229 (2012), arXiv: 1209.2332 [hep-ph] 

[64] T. Wolkanowski, Resonances and poles in the second Rie- 
mann sheet, Master’s thesis, Goethe-Universitat Frank¬ 
furt am Main (2013), arXiv:1303.4657 [hep-ph] 



